A sound card based multi-channel frequency measurement system 
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For physical processes which express themselves as a frequency, for example magnetic field mea- 
surements using optically-pumped alkali-vapor magnetometers, the precise extraction of the fre- 
quency from the noisy signal is a classical problem. We describe herein a frequency measurement 
system based on an inexpensive commercially available computer sound card coupled with a software 
single-tone estimator which reaches Cramer-Rao limited performance, a feature which commercial 
frequency counters often lack. Characterization of the system and examples of its successful appli- 
cation to magnetometry are presented. 



I. INTRODUCTION 

The present work is motivated by the need for a high 
resolution frequency measurement system for analyzing 
signals generated by optically-pumped cesium magne- 
tometers p] . A set of such magnetometers will be used 
for a detailed investigation of magnetic held fluctuations 
and gradients in an experiment searching for a neutron 
electric dipole moment (nEDM). The experiment calls 
for a magnetic field of between 1 to 2 /iT controlled at 
the 80 fT level when measured over 100 s time intervals, 
control corresponding to a relative uncertainty between 
40 to 80 ppb. The magnetometers are based on the fact 
that for low magnetic fields the Larmor precession fre- 
quency /l in a vapor of Cs atoms is proportional to the 
modulus of the magnetic field B 



Sl = 7 B 



(1) 



The proportionality factor 7 is a combination of fun- 
damental and material constants and has a value of 
« 3.5 kHz//iT for 133 Cs. The precession of the atoms 
modulates the resonant absorption coefficient of the ce- 
sium vapor, which is measured by a photodiode moni- 
toring the power of a laser beam traversing the atomic 
vapor . In the self-oscillating mode of operation 0, 
the magnetometer signal is of the form 



s(t) = Asm (2nf L t + <p) + s . 



(2) 



The Larmor frequency, has to be extracted from the 
signal. Equationnconnects the frequency determination 
precision directly to the resulting field measurement pre- 
cision. The basic demand on the frequency measurement 
system in order to achieve the required field precision is 
a resolution of a few hundred //Hz in an integration time 
of 100 s. Moreover, the synchronous detection of signals 
from an array of magnetometers requires a cost-effective 
multi-channel solution. 

In our recent study |3j of optically-pumped mag- 
netometer performance, frequency measurements were 
made with a commercial frequency counter (Stanford Re- 
search Systems, model SR620), which has a limited fre- 
quency resolution, thereby limiting the magnetic field de- 
termination. Frequency counters rely on the detection of 



zero crossings of a periodic signal in a given dwell time. 
Their performance is limited by their resolution of the 
zero crossing times, an event which is affected by the 
amplitude, offset, and phase noise of the signal. In de- 
manding applications, such as the one investigated here, 
that timing jitter limits the ultimate frequency resolution 
of the magnetometer signal measurement. Put simple, 
the limitation of frequency counters is due to the fact 
that they use only information in the vicinity of the zero 
crossings, while valuable wave form information from in 
between the zero crossings is ignored. 

As a more powerful alternative one can use numerical 
frequency estimation algorithms to extract the frequency 
from the complete waveform sampled at an appropriate 
rate and with a sufficient resolution. The performance of 
an ADC-based measurement system for measuring a sin- 
gle frequency of about 8 Hz was discussed in Q . Under 
the assumption that a stable clock triggers the ADC, the 
authors in 4| show that the lower limit of the frequency 
resolution of their system coincides with the Cramer-Rao 
lower bound (CRLB) @. The CRLB is a well-known 
concept from information theory and describes principle 
limits for the estimation of parameters from sampled sig- 
nals. 

In our application, the Larmor frequency in a mag- 
netic field of 2 fiT lies in the audio frequency range 
[Jl = 7 kHz). We have investigated whether a com- 
mercially available (and rather inexpensive) professional 
multi-channel sound card would present a viable solu- 
tion for sampling the magnetometer signals. The estima- 
tion of the frequency from the sampled data was done 
by a software algorithm. In the following we will show 
that such a simple system can indeed be used for CRLB 
limited real-time frequency measurements and for a de- 
tailed study of noise processes which limit the precision 
of atomic magnetometers. 



II. THE SYSTEM 

The frequency measurement system consists of a pro- 
fessional sound card (M-Audio Delta 1010) for digitizing 
the analog input data, an atomic clock to provide a stable 



2 



time reference, and a standard personal computer (PC) 
which reads the data and runs the frequency estimation 
algorithm. The sound card provides 8 analog input chan- 
nels in a breakout box that connects to a PCI interface 
card in the PC. The analog input signals can be sampled 
with a resolution of up to 24 bit at a sampling rate of up 
to 96 kHz. In order to limit the amount of data we used 
only 16-bit resolution, which was proven to yield suffi- 
cient precision. Jitter or drifts of the sampling rate in- 
duce additional phase noise on the sampled signal which 
can seriously degrade the precision of the frequency esti- 
mation. An essential feature of the Delta 1010 sound card 
is its "world clock" input which can be used to phase-lock 
the internal clock of the sound card to an external 96 kHz 
time base. The time base was realized by a frequency gen- 
erator synchronized to the 10 MHz signal of a rubidium 
frequency standard (Stanford Research Systems, model 
PRS10). The Rb frequency standard provides a relative 
stability of 10~ 12 in 100 s which minimizes possible sam- 
pling rate jitter and drifts far below the required level. 
The requirements for the PC system are not very de- 
manding as long as it allows for the continuous recording 
of the 16 bit data sampled at a rate of 96 kHz (5.8 GB/h 
for 8 channels). A 1.8 GHz Pentium-4 processor was fast 
enough for real-time frequency determination for all eight 
channels at a given integration time. However, for the de- 
tailed analysis described below, in which the integration 
time is varied, the time series were evaluated off-line from 
the stored sampled data. 



III. PERFORMANCE 

Considering the magnetometer signal given by Eq. 
the frequency Jl is to be determined from the AC- 
coupled signal data which, after sampling, are of the form 



x n =Asin 



k=l 



+ Sx n , (3) 



n=(0,...,JV-l), 



where A is the signal amplitude, At the time resolution 
(inverse of the sampling rate r s ), and 0o the initial phase. 
The number of sample points is N = t / At{= rr s ), where 
r is the measurement integration time. Also shown is 
the noise contribution at each point n arising from phase 
noise #</>„, frequency noise £^fc =1 <5/fe, an d offset noise 
5x n . The frequency is determined from the data x n by 
a maximum likelihood estimator based on a numerical 
Fourier transformation which provides a CRLB limited 
value [||- The algorithm iteratively searches for the fre- 
quency / that maximizes the modulus of the Fourier sum 



MF(f) = 



2 , x n W n exp ( i n 



(4) 



Under ideal conditions (stable field and ideal electron- 
ics), the frequency and phase noise (Sfk and 8(j) n ) are 
not present in the signal (Eq. EJ) . The fundamental noise 
contribution is the photocurrent shot noise, which is pro- 
portional to the square root of the DC offset sq in Eq. 
The noise is converted, by a transimpedance amplifier, 
to voltage V pc and has a Gaussian amplitude distribu- 
tion with zero mean, corresponding to a white frequency 
spectrum that is characterized by its power spectral den- 
sity p\ (in V 2 /Hz). The signal-to-noise ratio (SNR) is 
defined as A 2 /(plf s ), where f s = (2At)" 1 = r s /2 is the 
sampling rate limited bandwidth, i.e., the Nyquist fre- 
quency or the highest frequency that can be detected 
unaliased. The CRLB of the frequency estimation from 
such an ideal magnetometer signal is given by the vari- 
ance 



CRLB 



TT 2 A 2 T 3 ' 



(5) 



In frequency metrology it is customary to represent 
frequency fluctuations in terms of the Allan standard de- 
viation a a — or its square a A , the Allan variance 010. 
One can show that for white noise a a coincides with the 
classical standard deviation 9]. A double lo garithmic 
plot of the dependence of a a on the integration time r is 
a valuable tool for assigning the origin of the noise pro- 
cesses that limit the performance of an oscillator (see for 
example 0,0)- As shown in Table [I] the variance <j\ 
depends both on integration time r and measurement 
bandwidth f c , which, for a measurement interval r, is 
given by f c = (2r) _1 . When that relation between band- 
width and integration time is inserted into the formulas 
given in the central column of Table |U 8] , one finds the 
typical r dependencies of the Allan standard deviation 
aA shown in the right-hand column. In the presence of 



noise source 


°1{t) 


oa(t) 


white offset 


3 2 1 

„2 A 2 P* r 3 




flicker frequency 


2 In 2-h 2 r l 


ocr° 


white frequency 


1 2 1 
- 2 -Pf- 


OCT- 1 / 2 


flicker phase 


3 u2 ln(27r/ c r) 
4tt 2 r 2 


-1 

OC T 


white phase 


3 2 fc 


OC T-3/ 2 



where W n is a windowing function. 



TABLE I: The central column shows the dependence of the 
Allan variance a\ on the integration time r and measurement 
bandwidth f c for the noise sources listed at the left Q . White 
noise sources a are characterized by their power spectral den- 
sity p 2 . The frequency dependent spectral density of flicker 
noise process a is modeled by p 2 (/) = /i 2 //. By assuming 
the relation f c = (2r) _1 we find the power laws which typify 
each noise type, shown in the right hand column. 
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FIG. 1: Allan standard deviation of the frequency of a syn- 
thesized sine wave affected by different noise processes. From 
top to bottom: flicker frequency noise (black dots), white 
frequency noise (open circles), flicker phase noise (black tri- 
angles), white phase noise (open triangles), white offset noise 
(black squares). 



several uncorrelated noise processes, a, the variance of 
the estimated frequency is given by 



- 2 = E^(/) 



(0) 



Note that for a magnetometer signal, the contribution 
from Eq. [S] will always be present in the sum. 

We first investigated whether our data analysis algo- 
rithm reproduces the theoretical r-dependencies shown 
in Table ||| For that purpose we generated time series 
(16 bit, 96 kHz) corresponding to Eq. [21 with only one of 
the phase, frequency, or offset noise terms enabled, and 
selected with well defined spectral characteristics (flicker 
or white). Figure shows the Allan standard deviation 
&A of those synthetic data. The emphasis here lies on 
the slopes rather than on the absolute values, which were 
chosen to yield a readable graph. 

Next, we investigated the ability of the sound card to 
reach CRLB limited detection of a 7 kHz sine wave. The 
wave was generated by a digital function generator (Agi- 
lent, model 33220A) stabilized to the same Rb frequency 
standard as the sound card. In order to simulate a signal 
comparable to that of the magnetometers, the SNR of 
the function generator output was artificially decreased 
from its nominal value of better than 5 x 10 5 to about 
1.3 x 10 5 (in a 1 Hz bandwidth) by adding white offset 
noise. We recorded a 1 h time series of that signal, sam- 
pled with 16-bit resolution. The data were analyzed with 
the same algorithm as above and yielded an Allan stan- 
dard deviation cta(t), shown as black dots in Fig. [2Ja) - 
The measurement agrees on an absolute scale with the 
CRLB calculated using Eq.[5]and the applied SNR. In ad- 
dition to the offset noise, a second noise source was used 
to apply 1/ / noise, in turn, to the frequency or to the 
phase modulation input of the function generator. The 



1a 1 
10 2 
10 3 

N 
I 

< 10 4 



10 5 
~6 

"in 10 6 
c 

< 10 7 
10 8 
10 




















•a- a. ^-1 


frequency counter 






r 


0" -. . 1 






- 


r 


T 


'"•a'; A -A 








r 


''x 

sound card 




3 ----o-o-o a ---< 


) 

c) ] 










»>>.. 










A '"-A b) 

'% 

A _ 










a) ' 


1 



0.01 



0.1 1 10 

Integration time x (s) 



100 



FIG. 2: Allan standard deviation of the frequency of a sine 
wave affected by different noise processes, measured with the 
sound card (a-c) and a frequency counter (d): a) white offset 
noise, b) white offset noise and flicker phase noise, c) white 
offset noise and flicker frequency noise. The dashed lines rep- 
resent the dependencies calculated on an absolute scale using 
the applied noise amplitudes, d) The same signal as in a) 
measured with a commercial frequency counter (Stanford Re- 
search model SR620) with a 300 Hz input bandpass filter. 
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FIG. 3: a) ASD a a of magnetic field fluctuations inside a 
multi-layer magnetic shield, b) Residual fluctuations of the 
stabilized magnetic field. The dashed lines indicate the CRLB 
(left) and an assumed white frequency noise limitation (right), 
c) Stability required for the proposed nEDM experiment. 



resulting a a of the measured data is shown in Figs. Eh) 
and c). Figure [2Ji) shows a a derived from the same sig- 
nal as Fig.[2Ji) but analyzed by the commercial frequency 
counter (Stanford Research Systems, model SR620) that 
was used in The three points shown correspond to 
the three possible integration times of the SR620. It can 
be clearly seen that the counter technique does not allow 
the correct measurement of these faint noise processes. 
However, extrapolation of the data points suggest that 
for integration times less than 10 ms the CRLB could be 
reached. 

Finally, after the frequency estimator algorithm and 
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the sound card had proven their CRLB performance 
limit, we used the system to analyze the frequency gen- 
erated by an optically-pumped magnetometer (OPM). 
A magnetic field of 2 /iT was produced by a solenoid 
driven by an ultra-stable current source. The OPM sig- 
nal in that field is a 7 kHz sine wave. The OPM and the 
solenoid were located in a 6-layer magnetic shield in order 
to suppress external field fluctuations. Figure^) shows 
(7a of a 2 h time series recorded with the sound card. 
The data represent pure magnetic field fluctuations. In 
particular, the approximately 2 mHz fluctuations in the 
range between 2 to 200 s could be traced back to irreg- 
ular current fluctuations in the solenoid. Nevertheless, 
the relative field stability — and therefore the relative 
current stability — is on the order of 3 x 10~ 7 for that 
range of integration times. However, for a 100 s integra- 
tion time the field instability exceeds the requirement for 
the nEDM experiment mentioned in the introduction. 

In order to determine the magnetometer performance 
limit, we actively stabilized the magnetic field in the fol- 
lowing way. The magnetometer frequency was compared 
to a stable reference oscillator (i.e., the Rb frequency 
standard) by means of a phase comparator, and the er- 
ror signal was used to control the solenoid current, thus 
realizing a phase-locked loop. Figure 0b) shows the Al- 
lan standard deviation of the OPM in the stabilized field, 
which is CRLB-limited up to an integration time of 1 s. 
The noise excess between 1 and 300 s above the limits 
expected from the CRLB and the assumed white noise 



limitation shows the limitation of the current stabiliza- 
tion scheme, which nonetheless allows the suppression of 
the fluctuations by three orders of magnitude at the in- 
tegration time of interest. 

We have realized a frequency measurement system 
based on a digital sound card and have shown that it 
yields a performance superior to commercial frequency 
counters. We have proven that the system yields CRLB 
limited frequency resolution in measurements of sine 
waves affected by various sources of noise. We have used 
the system to prove that, at least in a limited range 
of integration times, an active field stabilization by an 
optically pumped magnetometer is limited by the the- 
oretical Cramer-Rao bound. The performance and the 
multi-channel feature of the sound card and its external 
frequency reference option present a low-cost alternative 
for applications requiring simultaneous characterization 
of several frequency generation systems, especially for 
long integration times. 
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